Lecture-1 D 
More Satellite/Sensor 
Geometry (Rev-C) 


EEL 5432 Sat Remote Sensing 



Lesson Objective 

To continue our discussion about satellite/sensor 
geometry as it applies to remote sensing 

To introduce the Earth’s latitude/longitude grid and 
“great circle” navigation 

Present orbit plots from Systems Tool Kit 



Great Circle Navigation 

Consider a spherical earth with a center “C” and two surface 
locations (“A” & “B”) 

Pass a plane through these three points, and the intersection with 
the sphere forms a “great circle” 

The shortest distance (between these points along the surface) is 
the arc AB 



The Earth’s Navigation Grid 

(Latitude & Longitude) 



Longitude lines are great circles formed by planes that 


contain the earth’s axis of rotation 


Longitude 

Lines of longitude 
run across the globe 
from North to South, 
but measure East and 
West starting at the 
Prime Meridian. 
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Latitude lines are formed by planes that are 
perpendicular to the earth’s axis of rotation 


Latitude 



Equator 0° 
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Lines of latitude go 
across the earth from 
East to West, but they 
measure the globe from 
North to South starting 
at the Equator. 
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Earth Navigation Grid: Latitude & Longitude 
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In this class, we will use System 
Tool Kit (STK) Software to make 
2D & 3D images of Low Earth Orbit 
(LEO) Satellite Ground Tracks 



Review of Orbit Elements 



Orbital Elements 
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Orbital Elements 



View AA - side view 
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Tropical Rainfall Measuring 

Mission (TRMM) 


Altitude: 399 km 
Inclination: 35.0 deg 

(prograde orbit = inclination ^ 90°) 

Orbit period: 92.5 min 



TRMM Sub-satellite point “ground track” 
for single orbit (period = 92.5 minutes) 

Rule of thumb” all LEO orbits have periods ~1.5 hours 









TRMM orbit ascending node 
- the point where the satellite passes through 
the equatorial plane heading North 












TRMM - satellite direction of flight at the 
ascending node is equal to the orbit inclination 














TRMM ground-track max (minimum) latitude 

= the orbit inclination 












TRMM orbit descending node 
- the point where the satellite passes through 
the equatorial plane heading South 


















TRMM - satellite direction of flight at the 
descending node = negative orbit inclination 












rRMM - Note that the ascending nodes for 
consecutive orbits are displaced westward 

why? 










Example: Spherical View of North Pole 
@ time = T 0 start time (asc node) & Sat Lng = 0 


Earth’s Rotation 



Lng = 0° 







Example: Spherical View Viewing North Pole 

During one orbital period (92.5 min) the earth rotates CCW, 
so the starting longitude of the next orbit moves westward 


Earth’s Rotation 












TRMM - single orbit (92.5 minutes) 








TRMM - 2 orbits 












TRMM - 7.8 orbit’s (12 hours) 












TRMM - 15.6 orbit’s (24 hours) 





QuikSCAT 

Altitude: 804 km 
Inclination: 98.6 deg 

(retrograde orbit = inclination > 90°) 

Orbit period: 100.9 min 




QuikSCAT - single orbit (101 min period) 



























QuikSCAT - Asc & Desc Orbit Segments 
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QuikSCAT - 3 orbits 
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QuikSCAT - 7 + orbits (12 hours) 














QuikSCAT - 14.3 orbits (24 hrs) 














TRMM & QuikSCAT 
Collocations 



QuikSCAT & TRMM Satellite - 1.5 Orbits 
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1 st QuikSCAT & TRMM Collocation 
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2 nd QuikSCAT & TRMM Collocation 
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3 rd QuikSCAT & TRMM Collocation 



















4 th QuikSCAT & TRMM Collocation 
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Summary: Orbital Properties 

The maximum latitudinal extent of the 
subsatellite point is “/ “ (orbit inclination) 

- If / < 90°, the orbit is prograde, and the rotation of 
the satellite, projected on Earth’s equatorial plane, 
is in the same direction as the rotation of the Earth 
(eastward) 

- If / > 90°, the rotation is retrograde (opposite 
rotation = westward) 

- If / = 90°, the satellite is in polar orbit; and 

- If / = 0°, the orbit is equatorial 

The ascending node N is the intersection of the 
subsatellite track with the equator at the time the 
satellite crosses the equatorial plane going northward 



Sensor Spin Rate Calculation 



Conical Scan: Measurement Arc (i^) 



U Center of 
the Earth 






Conical Scanning Sensor Geometry - 
Measurements Contiguous Spatial Sampling 




Contiguous Spatial Sampling “just touching” along the scan 
direction and scan-to-scan in the direction of flight (along the 
satellite sub-track) 




Conical Scanning Sensor Geometry - 
Measurements Nyquist Spatial Sampling 



Nyquist Spatial Sampling 50% overlap along the satellite sub¬ 
track and contiguous sampling along scan 




Conical Scanning Sensor Geometry - 
Measurements Spatial Sampling - cont.-l 


Contiguous Spatial Sampling 
along track & Nyquist spatial 
sampling cross track (50% 
overlap) 



Nyquist Spatial Sampling 
along & cross-track (50% 
overlap) 





Conical Scanning Sensor Geometry 



T 

Sat. Velocity 

Polar View 






















Conical Scanning Sensor Geometry cont. 


Cross-Track Plane 

sin(y ct ) = sin(y) * sin( W) 
y ct = sirr^sir^y) * sin( *^0} 

where Y is measured in incident plane 

Swath Width =2 p y ci 


Swath Width 



Note: y ct is in units of radians 


Note: Swath width is measured in the cross-track plane 











Conical Scanning Sensor Geometry cont-2 

Swath Width = 2 py a Note: y ct is in units of radians 

where: 

p = earth's radius = 6378.155 km 

y a = central angle of swath in the cross - track plane 

= arcSin (Sin ( Q - <E> )* Sin (ip ) ) 

jp = half azimuth of the active (measurement) swath 
Q x = incidence angle 
O = nadir (cone) angle 

Y = ( 0 ~ °) 



Conical Scanning Sensor Geometry cont-3 
(Contiguous beams along-scan) 



Active scan length, S = 2 * ip * r 
r = R* sinO and ip is in radians 


# Beams = Int 
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Example: Sensor Calculations Conical imager 

(Incident Plane) 

h = 405.0 km, O = 48.0°, p = 0.624° = (0.01089 rad), W = ± 65.0° 

Incidence angle 
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y =6. — d) = 4.22° (0.0736rad) 
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Slant range 


R = p 


( • \ 

sin y 

ysin0 ) 


= 631.11 km 


IFOV along subtrack: 

along - track = ^ &/ = 11.22 km 

/cos 6 

cross - track = R - (3 = 6.87 km 



Example: Conical imager Calculations cont. 


Active scan length, S = 2 * ip * r 
r = R* sinO and ip is in radians 

S = (2 x 65° x pi/180) x 631.11 xsin(48°)= 1064.18 km 


# Beams 


Int 


IFOV 


cross-trk 


+ 1= 155 


# of contiguous IFOV’s in the active arc length 



Example: Conical imager Calculations cont.-2 

Calculate antenna spin rate for contiguous sampling (along-trk) 



(3.986^ + 5 
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1.2\km / s 


During one scan, the along-track distance moved is 

Ax = IFOV along _ track = 11.22 km 


The required spin period is 

T s = Ax/V„ = 11.22/7.21 = 1.556 sec 

The antenna spin rate in revolutions/minute is 
Spin rate = 60/T s = 60/1.556 = 38.55 rpm 



Example: Conical imager Calculations cont.-3 

Calculate antenna spin rate for Nyquist sampling (along-trk) 



(3.986^ + 5 
V 6783.16 


7.67 km / s 


V 

g 



_ „ 6378.155 
7.50*- 

6783.155 


1.2\km / s 


During one scan, the along-track distance moved is 

At = (IFO V along _ track )/2 =11.22/2 = 5.61 km 


The required spin period is 

T s = Ax/V„ = 5.61/7.21 = 0.778 sec 

The antenna spin rate in revolutions/minute is 
Spin rate = 60/T s = 60/0.778 =77.11 rpm 
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Calculate Swath Width Req’d to fill the “Gap” at the Equator 

(worst case) 


STK View 1 - Earth 
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Swath cover gap 










































Calculate Swath Width for Contiguous Imaging at 

the Equator (no swath gap) 
































































































Example: Swath required to “fill the gap” @ equator 


Consider a satellite with the following parameters: 

Circular orbit h = 405 km ; orbit inclination i = 65° 

Find the sensor swath to provide 100% earth coverage? 

Solution: 

Orbit period from Kepler’s 2 nd law = 92.6633 min = 5559.8 sec 
Earth’s angular rotation rate = (2 pi)/sidereal day= 2 pi /86164.0905 

= 7.292xl0 -5 rad/sec 

Earth’s rotation/orbit @ equator = Orbit_period*earth_rotation_rate 

= 0.40527 rads = 23.22 deg 

Orbit sub-track separation = (earth’s rotation/orbit) * earth’s radius 

= 0.40527*6378.155 = 2584.84 km 
SW = separation*sin(i) = 2582.27 x sin(65°) = 2342.66 km 



Example: Swath required to “fill the gap” cont. 


Consider a satellite with the following parameters: 

Circular orbit h = 400 km ; orbit inclination i = 65° 

Req’d SW = 2342.66 km 

Assuming max Az angle of ± 90°, what is the incidence angle 
required? 

Solution: 

For this case (Az = 90 deg), the SW is like a cross-track scanner, so 
we can solve for the earth central angle (gamma) that produces the 
required S W = 2 *y*p or 

y = SW/2p radians = 2342.66/(2*6378.155) 

= 0.18365 radians (10.52 deg) 



Example: Swath required to “fill the gap” cont. 


Consider a satellite with the following parameters: 

Circular orbit h = 400 km ; orbit inclination i = 65° 

Req’d y = SW/2p radians = 0.18365 radians (10.52 deg) 

R p + h 

From the law of sines t = . 7777 

sin y sm(fcf) 

From the law of cosines 

R = \p 2 + (p + hf -2 p(p + h ) cosyl 

Solve these equations simultaneously for R and 6 yields 
R = 1272.45 km & 6= 76.78° & 0= 66.26° 

Note that the max cone ang for the line of sight (LOS) to be tangent 
to the limb of the earth: <p max = Arc Sin (p/(p+h)) = 70.10°; thus the 
above solution is valid 



